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Abstract 

The universal curve a /A of nuclear photoabsorption is investigated 
within a Fermi gas model of nuclear matter. An energy range from pion 
threshold up to 400 MeV is considered. The interactions between nu- 
cleon, pion, A-isobar and photon are considered in the non-relativistic 
approximation with corrections of the order 1/M taken into account 
with respect to proton mass. Analytical expressions are obtained, in 
which the influence of nuclear correlations, two-nucleon contributions 
and relativistic corrections is studied explicitely. An extension of the 
model calculation to nucleon knock-out reactions is discussed. 
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1 Introduction 



The phenomenology of nuclear photoabsorption is governed by two char- 
acteristic features. First, all nuclei with mass numbers A ranging from 10 
to more than 200 obey the same fundamental curve a(u)/A for the total 
photoabsorption cross section devided by A as a function of the photon en- 
ergy oj. Second, the A-isobar excitation of the nucleon is responsible for 
the main properties of this curve in the energy region between 200 and 
400 MeV. Models, which focus on the behaviour of the A-isobar in a nu- 
clear environment, namely the A-hole calculations (2|, have proven to be 
highly successful in explaining the experimental findings for pion scattering 
processes ||. Indeed, on grounds of the A-hole formalism a wide variety 
of pion-nucleus reactions can be described within one consistent framework 
||. In the case of photonuclear reactions, however, serious descrepancies 
remain, which partially have been accounted for by including non-resonant 
background terms H, p|. Nevertheless, such a procedure, in particular for 
nuclear photoabsorption, either leads to contradictions with previous A- 
hole results or lacks the complete agreement with experimental data ||. 
The question arises, whether theoretical ingredients other than in- medium 
A-hole propagations can lead to a similar degree of accuracy. Therefore 
it is natural to address the situation from a different point of view asking 
to what extend the absorption process can be described, when only a very 
simple A-nucleon interaction is used and all additional effects are accounted 
for in a purely diagrammatic approach. When combined with a simple form 
of nucleon momentum distribution inside the nucleus, namely a Fermi gas 
model, this leads to analytical expressions, in which different corrections 
to this lowest-order approximation can be studied. This is the aim of the 
present article. In our formalism we follow closely Wakamatsu and Mat- 
sumoto @. However, we do not introduce a phenomenological potential to 
account for the deviation of the nucleon wave functions from plane waves, 
but study the influence of such corrections in a perturbative way, similar to 
our previous work j8|. In addition, our focus is on the energy-dependence 
of the total photoabsorption, rather than on the differential cross section 
as a function of the momentum of the outgoing proton. A characteristic 
feature of Wakamatsu's and Matsumoto's approach is the technically equal 
treatment of the (7,pn) and the (7,pp) knock-out process, which allowed 
them to obtain a natural explanation for the supression of the two-proton 
knock-out. This characteristic property is also present in our calculation, 
where it is related to a vanishing trace in spin space. In the last years 
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the ratio (7,pp) / (7,pn) has been investigated in depth within different for- 
malisms. In an extension of Wakamatsu's and Matsumoto's work by Boato 
and Giannini ||, finite-size effects have been calculated and, more recently, 
a combination of pion exchange and shell-model wave functions was used to 
investigate this quantity |10| ]. Currently, two complementary approaches for 
the description of nuclear knock-out reactions exist. Carrasco and Oset jy]] 
used a diagram-oriented many-body expansion in a Fermi gas. The evalu- 
ation of self-energy diagrams leads to an accuracy for medium effects high 
enough to study knock-out reactions in great detail. With their primary 
goal being thus different from ours, their formalism does not yield isolated 
expressions for the resonant and non-resonant parts of the mechanisms of 
nuclear photoabsorption A different approach is used by the Gent group 



[12, [KJ, |i4| , where the main emphasis lies in the construction of realistic 
shell-model wave functions, rather than on a microscopic description fully 
based on the evaluation of Feynman diagrams. As the nuclear photoabsorp- 
tion is almost insensitive to structural differences between nuclei, the quality 
of their approach becomes obvious in the investigation of differential cross 
sections for nucleon knock-out, rather than of photoabsorption. 

Nuclear photoabsorption provides an interesting tool to study the interplay 
between one-nucleon and two-nucleon contributions. We obtained analyt- 
ical expressions for these contributions, as well as for their resonant and 
non-resonant parts. This set of results can be used as a starting point to 
include (and test) further nuclear or nucleonic effects. In Section 2 the basic 
notations are listed, as well as the interaction terms and the most important 
model assumptions and approximations, which are present in this calcula- 
tion. The main results and their most prominent properties, e.g. the effects 
of relativistic corrections and nuclear structure, are discussed in Section 3, 
where also the following possible extension of such an approach is considered: 
As the angular dependence of the two-nucleon process is not very strong (cf. 
p5[), the different mechanisms contributing to the photoabsorption curve 
can also be used to understand qualitative features of experimental data for 
nucleon knock-out processes as a function of the photon energy. In Section 4 
some concluding remarks are made, with an emphasis on the applicability of 
the partial cross sections, whose analytical forms are given in the Appendix. 
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2 Formalism and Notation 



The starting point of our investigation is the static Hamiltonian for the 
pion-nucleon interaction, 

/ 

H-kNN = a-Vr-cf) (1) 

m — 

together with a minimal coupling to the photon field. In eq. dlj) m is the 
pion mass. For all coupling constants we use the same notation and values 
as given in [16|, in particular / 2 /(4-7r)=0.08. In eq. (|l]) underlined symbols 



denote vectors in (cartesian) isospin space, while an arrow indicates a vector 
in coordinate space. In the static limit the interactions with the A-isobar 
excitation of the nucleon are determined by the following Hamiltonians (see 
e.g. @): 

H = S + .(Vx A) T+ (2) 

and 

H nNA = ~— (3) 
m — 

with the hermitian conjugate to be added in both cases. Here S and T are the 
l/2-to-3/2 transition operators in spin space and isospin space, respectively; 
e is the proton charge, e 2 = 1/137. For the coupling constants we have /a=2 
and / 7 7va = 0.35. In all cases, where high momentum transfers occur at 
the pion-nucleon vertices we regularize the vertex functions by introducing 
dipole form factors 

. A 2 -m 2 

M?) = A 2_ q 2 ( 4 ) 

as was also done e.g. in 0, 17\. The value for the cut-off parameter A 



has been taken to be 800 MeV. This value gives the best agreement of our 
predictions with the experimental data. In addition, a similar value has 
been used in |L8], [7j]. The general expression for the total cross section (J\[u) 
of the photoabsorption with one nucleon outside the Fermi sphere and one 
pion in the final state (one-nucleon process) is of the form 

d< l 4?F | Tl |2 

2e (? (2vr) 3 2uj 1 1 




- 1 n(p) 1 — n(k + p — q) . (5) 



(k+p 



2M 
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For the total cross section 02 (^-0 of the process with two free nucleons in the 
final state it is given by 

r VdpiV dp 2 r dp 3 dp4 47T 2 
J (2vr)6 J (2vr)6 2u 1 2| X 

n(pi)n(p 2 ) [1 - n(p 3 ) ] [1 - n(p 4 ) ] . (6) 

The notation for the external momenta is shown in Fig. ((!]). The function 
n(p) = 9{pf — \p\) is the occupation number, 9{x) is the step function. 
Furthermore, V is the nuclear volume, pf is the Fermi momentum, M is 
the mass of the proton, V is the nuclear volume, V = 3ir 2 A/(2p F ), and 

£ q = \J\q\ 2 + m 2 is the energy of the outgoing pion. 

Both amplitudes T\ and T2 consist of non-resonant and resonant parts, Tj = 
t!> NR ^ + T± . Diagrammatically this decomposition is shown in Fig. (^). 
The second (crossed) contribution to the resonant part is small due to the 
big energy denominator and will be skipped in the following. In the energy 
5-function in eq. (|) and eq. (H) we will usually not consider the smallest term 
connected with the kinetic energy of the incoming nucleon, as it is smaller 
than p 2 F /2M » 37MeV. Indeed, as we expect its contribution to introduce 
only a small modification of the actual p-dependence in the integrand in 
(||), we substitute it by its average value (p 2 ) /2M = (3/5)p 2 F /2M, which 
results in an overall shift of the absorption cross section. For the one- 
nucleon process we find first-order relativistic corrections to be essential for 
a successful treatment of the absorption process. In the case of the resonant 
contribution, such corrections are accounted for by making the following 
substitutions in the vertices JjJ: 

~~ (p + k), (7) 



k — > k(l + — ) - — p, (8) 
\ MJ M v ; 

where is the A-isobar mass and A = — M. For the non-resonant 
part, the corrections give amplitude T\ of the following form: 



rp(NR) V2ef . . . 

m ' 



1 + -^- 
2M 



(q — k) 2 + m 2 — (e q — uj) 

2i(a-q)(e-(k+p-q)) 



(k+p- q) 2 -2Mlo - (p 



2i(a-(q-k))(e-q) 

-to) 2 

0) 
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which corresponds to the production of a ir~ . In comparison with j^] we 
neglected those terms in (g), which modify the result by less than 2 per cent. 
It should be noted that no free parameters are present in our approach. A 
certain model dependence exists, however, in the selection of diagrams. We 
neglect all those diagrams, which are suppressed by some mechanism. In 
Fig. (|3D two examples for suppression mechanisms are given. For Fig. (|a) 
the contribution is small because in the case of infinite nuclear matter due 
to momentum conservation the four-momentum of the photon should be 
equal to that of the outgoing pion, which is impossible. For Fig. (||b) let us 
consider the case, where the upper two nucleons (incoming and outgoing) 
are identified. Furthermore, let us select a A-isobar as an intermediate state. 
Then, one has a vanishing trace in spin space: 

Sp \(S + ■ (k X £)) {S-q}] =0 

at q = k. Therefore, the contribution in this case vanishes in the non- 
relativistic limit considered here. As a result of applying such methods to 
the various diagrams, we have obtained that only those displayed in Fig. (|l|) 
should be taken into acount. 



3 Results 



Explicit evaluation of the diagrams shown in Fig. (Q) leads to amplitudes for 
the one- and two-nucleon contribution to nuclear photoabsorption. Squaring 
these amplitudes and summing over spin and isospin states of the nucleons 
by using trace methods as previously || one finds the expressions |Ti| 2 and 
| T2 1 2 , which enter into eqs. (||) and (^). It has turned out to be convenient 
to investigate the resonant and non-resonant parts of each of these contri- 
butions separately. This can be done by neglecting the interference terms, 
which are highly suppressed (cf. Fig. (0)). Shown here exemplary for the res- 
onant parts, one obtains the following expressions, which serve as a starting 
point for the integrations with respect to nucleon momenta: 

4e 2 // A / 7 iVA x2 
27iru V 



a 



(«), 



Vdp 



J dq x 



9 w r 
" + M 



(p~ kf -(p- q)' 



m 



n(p — k) [1 — n(p — q)\ 



(QxK)-e +Q 2 (Kxe)' 



(lo- A-p 2 /2My + F 2 /A 



(10) 
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where 



and 



uj -> ( A\ A 

Q = «-M~J> K = k { 1+ M)-M P 



a y 2 >{uj) 



16e 2 / /a/ 7 ata 



8l7r \ m 
V dpi V dp2 



UJ 



(27T) 
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(w - A) 2 + r 2 /4 
3 dpi n(pi)n(p 2 ) [1 - n(p 3 ) ] [1 - n(p 4 ) ] x 



2a 2 [a 2 + 3(e-a) 2 ] 4 , 2w 2 (e • a)(e- 6) 2 2 1 



(a 2 + m 2 ) 2 



(a 2 + m 2 ){b 2 + m 2 ) 



with 



a = P4 - P2 



Pi -P3- 



The analytical expressions for the resulting partial cross sections are given 
in the Appendix. For the total absorption cross section cr(u) a compari- 
son with experimental data is shown in Fig. ([|). The result of this model 
calculation compares favorably with the data. The one-nucleon and two- 
nucleon parts of the cross section are equally important at energies around 
250 MeV. As can be seen in this plot, significant features of the data, e.g. 
the position of the peak, are only obtained due to the interplay between the 
two mechanisms. As mentioned before, each of the contributions to the full 
curve in Fig. (|j) has a resonant and a non-resonant part. In Fig. (||) and (||) 
this decomposition is shown for the one-nucleon and the two-nucleon mech- 
anism, respectively. Here it is clearly seen that the non-resonant parts give 
an important contribution at lower energies. In the two-nucleon case the 
non-resonant part decreases with energy, while in the one-nucleon process 
it remains almost constant. 

It is interesting to see, in what way the one-nucleon partial cross sections 
are affected by the use of only the static limit of the interaction. Neglecting 
the first-order relativistic corrections in the current one has 



a 



~(^) 



Ae 2 (f A f lNA )< 



Qmax 



27m 2 

2' 



qdq (1 + S (g))x 



m 



~M 



+ 



1+ M ) +m 




(12) 
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and 



~af R \u) = ^L ( qdq(l + So(q)) 



4q 2 uJ 2 



2m?q 2 

v 

(q 2 + m 2 ) 2 

(13) 



where the integration limits in both cases are given by 



i mm, max 



MJ 

and the integrand contains the function 




In Fig. (||) the corresponding cross sections are compared with those resulting 
from eq. ( |To| ) and its non-resonant counterpart. The relativistic corrections 
lead mainly to a shift of the one-nucleon curve. This effect is essential 
for obtaining a good agreement with the experimental data. Although no 
relativistic corrections in the current are included, note that in ( |l2| ) and 
( |l3|) terms of that order have been kept in the kinematical contributions to 
the integrand. As the two-nucleon mechanism gives a comparatively small 
contribution at higher energies, we neglect relativistic corrections in 02(0;). 
This has also been done in 

A more difficult problem is the influence, which nucleon correlations inside 
the nucleus can have on the nuclear photoabsorption process. We investigate 



this aspect for the non-relativistic forms (12) and ( |13D of the one-nucleon 
case. The main reason for doing so is the fact that due to the square of the 
amplitudes we can express part of the integrand via the standard lowest- 
order central correlation function of a Fermi gas (see e.g. JOJ). The object, 
which is obtained by diagrammatically squaring the one-nucleon contribu- 
tion of Fig. (||a), can be coupled to an additional nucleon. In the incoherent 
case of the diagrammatical square also a further pion exchange can be al- 
lowed. The modifications of the correlation function, which occur due to 
such effects, have been investigated analytically in []|], where a corrected 
central correlator Sc(q) has been constructed. The function Sc is given in 
Fig. (H|). By making the substitution So(q) — ► Sc(q) these further corre- 
lations can be incorporated effectively. In Fig. ( |i~0| ) the one-nucleon contri- 
bution resulting from Sc is compared with the original form, in which So 
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has been used. It can be seen that such medium effects modify the result 
by about 15 per cent. Naturally, the use of this substitution technique can 
only be used to obtain an estimate for such a medium-induced modification. 
A full examination should involve the inclusion of the additional two-and 
three-nucleon diagrams in eq. (S). 

The limit of our approach is certainly reached, when a comparison with 
differential cross sections is attempted. An extreme case is the comparison 
with 4 He, for which a recent measurement of both, the one-nucleon and the 
two-nucleon channel exists |2C]. We compare the calculated average cross 
section with the data for the differential cross section in cm. frame, as 
we expect the angular dependence not to be strong in that frame. The 
corresponding plot is shown in Fig. (|iT|). Although no full agreement is 
obtained, it is interesting to note that the general features of the two cross 
sections are well reproduced, such as the peak positions and the relative 
size of the two processes. By these means it is possible to unambiguously 
identify the physical mechanisms behind the data points. A similar degree 
of agreement is obtained for other data plL 



4 Conclusion 



In the present paper we have developed a diagrammatical description of 
the nuclear photoabsorption process. The main result of our investigation 
is that the total photoabsorption cross section can be fully understood in 
terms of a simple physical picture, where point-like nucleons and A-isobars 
interacting via pion exchange are the relevant degrees of freedom. Due to 
the diagram-oriented formalism and the Fermi gas model as an approximate 
description of the nucleons in momentum space, we could obtain analytical 
expressions for all the relevant contributions to the photoabsorption curve. 
In this way a flexible and efficient description has been obtained, which 
can be used as a starting point for the investigation of additional effects. 
Especially in the low-energy part of our calculated curve, the agreement 
with experiment comes about as a non-trivial interplay between the one- 
nucleon and two-nucleon contributions. It is worth noting that, as long as a 
comparatively low cut-off parameter in the vertex form factor is used, there 
seems to be no need for an explicit diagrammatical inclusion of the p-meson 
as an additional mechanism of the nucleon-nucleon interaction. We found 
relativistic corrections in the case of the one-nucleon process to be crucial for 
obtaining a good agreement. The aspect of additional nucleon correlations, 



S 



which can be accounted for as a deviation of the nucleon wave functions 
from plane waves, deserves some further attention in future investigations. 
We could estimate the overall effect to be of the order of 15 per cent. 
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Appendix: Analytical expressions for absorption 
cross sections 



Here we present the explicit expressions for the four contributions to the 
photoabsorption curve, which have been used to obtain the figures shown 
in Section 3. As was mentioned earlier, the interference terms between 
the resonant and the non-resonant contributions are small (cf. Fig. (|7|)). 
Therefore, we can write the absorption cross section as a sum of four parts, 

ct(cj) = a[ NR \(v)+a[ R \u) + a2 NR \io) + a2 R \uj)- First, we deal with the one- 
nucleon case. In the case of the non-resonant contribution the absorption 
cross section can be represented in the following form: 



a 



(NR) 



3Ae 2 



w+ q+PF 



Am 2 <jj 2 Pp 



dq / dpgl(q)p[p 2 F - (q - p) 2 ] x 



G(p, q, oj) 6{2qoo — m 2 — q 2 



M 



P 



where 
G(p,q,u) 



1+— (— + —* 
2MJ D\D 2M 



(16) 



(17) 



and l{q) = Max[p^, \q— pp\] , D = q 2 +m 2 , uo± = uj±^/lo 2 — m 2 . In eq. (li~6|) 
the form factor g 2 (q) is the same as in (|j). Note that the integration with 
respect to the variable p can easily be performed, but the result is too lengthy 
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to be given here explicitely. The resonant part of the one-nucleon process 
has the following form: 



-1 



(pf ~ P 2 (1 + oj/M) 2 -u 2 + 2ujpy (1 + u/M) 



H(p,y,u) (18) 



[uj - A - p 2 ( 1 + u/M) 2 /2M A ] 2 + T 2 /4 
where b = (uj + pp)/(l +u)/M) and the integrand is given by 

i?0, y, uj) = 2A l [9(F - p - p F ) + 0(F) 0(p - p F - F) ] + 
1 



(x + l)Ai - -(x - x J )A 2 ]0(F - |p - PF |) 0(p + PF - F) , (19) 



with 



F 



1 _ 2py \ 



m 



-,1/2 



, X 



p 2 + F 2 - p\ 
2pF 



M J 1 + u/M 

Ai=u> 2 - 2toay + a 2 , A 2 = lo 2 (1- 3y 2 ) + 4way - 2a 2 
_ pA n+co/M\ 



(20) 



The mass difference A between the proton and the A-excitation is A = 292 
MeV, the width T of the A-isobar has been taken to be 115 MeV. Again, 
the integration with respect to y in eq. ( |l8| ) can be performed analytically, 
but due to its length the result is not presented here. For the partial cross 
sections of the two-nucleon case, we obtained the following result: 



o\ >(uj)+a 2 '(uj) 
L 2 (Q) 



(21) 



6(uj-E)6(5e-uj) 



L 3 (Q) 



$(/?!, Q)+ J $(02, Q) 



+ 



0(u - 5e) 6(9e -co) J $(/3i, Q) \gl(p) [g {nr) (p,uj) + G^ R \p,oj)] dp, 

Li(Q) ' 
where Q = VMuj , e = p 2 F /M , 

a P + PF a PF 

Pi = arccos — -= — , p2 = arccos 



V2Q 



V2Q 
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The elementary function $(/?, Q) is 

*(A Q) = U(p 2 "Pf)(2Q 2 +P 2 -^)sin 2 x + Q 4 (sin 4 x- ^sin 6 x) + 



1 



r l 



~ 2 Q 2 {x - - sin 4s ) - v^pQ [-(2Q 2 + p 2 - p%) sin 3 x - -Q 2 sin 5 x] + 



[^(2Q 2 + P 2 - P 2 F ) cos 3 x - Iq 2 cos 5 
o 5 



5 

x=tt/2-P 
x=/3 



(22) 



The functions in the integrand of eq. fl22|) , which characterize the resonant 
and non-resonant part, are given by 



p 



and 



8e'Afj NA M' 

27ir 2 m 6 p F (uj - A) 2 + T 2 /4 " (p 2 + m 2 ) 2 

2e 2 vlM 2 



7r 2 m 4p3 i 



1 + 



m 



(p 2 + m 2 ) 2 

respectively. In eq. ( f22D the integration limits are 



(p 2 + m 2 ) 2 



(23) 



(24) 



L 1 (Q) = Q-p J r, L 2 (Q) = J2Q 2 -p 2 F -p F 



and 



i 3 (Q) = V^ 2 -^ 



n 
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Figure Captions 



Figure 1 Notation for three-momenta of the external particles a) for the 
one-nucleon process and b) for the two-nucleon reaction. The wavy lines 
denote photons and dashed lines denote pions. A circle indicates a bound, 
an arrow a free nucleon. 

Figure 2 Diagrammatical forms of the resonant part (first two terms ) 
and the non-resonant part (last term) of the amplitude Tj. These 

contributions to the jit N N -interaction enter into the diagrams shown in 
Fig. g). 

Figure 3 Examples of diagrams, which have not been considered in this 
approach (cf. discussion in the text). 

Figure 4 Comparison of the calculated curve a(ui)/A for nuclear photoab- 
sorption with the experimental data. The dotted curve is the one-nucleon 
contribution o~\(uj), while the dashed curve represents the two-nucleon mech- 
anism o~\(oS). The data are taken from Ref. [2qJ . The empty (full) circles 
correspond to 208 Pb ( 12 C) data, while the squares represent data on 16 O. 

Figure 5 Resonant (dashed) and non-resonant (dotted) contribution to the 
one-nucleon reaction (full line). The corresponding analytical expressions 
a[ R \u>), a[ NR \u>) and o~\(uS), respectively, can be found in the Appendix. 

Figure 6 Same as Fig. ffi), but for the two-nucleon reaction. The analytical 
expressions (uj), (lo) and 02 (w), are also given in the Appendix. 

Figure 7 Contribution of interference terms. The full (dashed) curve cor- 
responds to the one- (two-) nucleon case. As can be seen from the overall 
scale, both for o\ and 02 this contribution is highly suppressed. 



14 



Figure 8 Effect of corrections of the order 1/M in the one-nucleon mech- 
anism. The full (dashed) curve is with (without) relativistic corrections. 



Figure 9 Correlation function Sc- The dashed curve corresponds to So, 
while in the dotted curve further two-nucleon correlations and three-nucleon 
correlations have been included. 



Figure 10 Effect of two- and three-nucleon correlation in the case of the 
non-relativistic one-nucleon part of the photoabsorption cross section. The 
dashed curve contains only So, while in the full curve Sc (cf. Fig. (\Q)) has 
been included. 



Figure 11 Approximate description of differential cross sections for he- 
lium. The differential cross section with respect to the direction of the outgo- 
ing proton is shown as a function of the energy of the incoming photon. Data 
points are taken from J1Z|/ . The full curve and filled circles correspond to 
the one-nucleon process, while the dashed curve and empty circles represent 
the two-nucleon case. 
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